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Abstract 

We show that doubly peaked electric fields are necessary to describe grazing-angle charge collection measurements 
of irradiated silicon pixel sensors. A model of irradiated silicon based upon two defect levels with opposite charge 
states and the trapping of charge carriers can be tuned to produce a good description of the measured charge 
collection profiles in the fluence range from 0.5x10^* Ucq/cm^ to 5.9x10^* Ueq/cm^. The model correctly predicts 
the variation in the profiles as the temperature is changed from — 10°C to — 25°C. The measured charge collection 
profiles are inconsistent with the linearly-varying electric fields predicted by the usual description based upon a 
uniform effective doping density. This observation calls into question the practice of using effective doping densities 
to characterize irradiated silicon. 
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1. Introduction 

A silicon pixel detector [1] is currently being 
developed for the CMS experiment at the CERN 
Large Hadron Collider (LHC). The detector will be 
a key component in the reconstruction of primary 
and secondary vertices in the particularly harsh 
LHC environment characterized by large track 
multiplicities and high radiation backgrounds. 
The innermost layer, located at only 4 cm from the 



beam line, is expected to be exposed to a 1 MeV 
neutron equivalent fluence of 3 x 10^* Uoq/cm^ per 
year at full luminosity. 

The response of the silicon sensors during the 
detector operation is of great concern. It is well un- 
derstood that the intra-diode electric fields in these 
detectors vary linearly in depth reaching a maxi- 
mum value at the p-n junction. The linear behavior 
is a consequence of a uniform space charge density, 
iVoff, caused by thermally ionized impurities in the 
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bulk material. It is well known that the detector 

characteristics are affected by radiation exposure, 
but it is generally assumed that the same picture 
is valid after irradiation. In fact, it is common to 
characterize the effects of irradiation in terms of a 
varying effective uniform charge density. In [2] we 
have proved that this picture does not provide a 
good description of irradiated silicon pixel sensors. 
In addition, it was shown that it is possible to ade- 
quately describe the charge collection characteris- 
tics of a heavily irradiated silicon detector in terms 
of a tuned double junction model which produces a 
doubly peaked electric field profile across the sen- 
sor. The modeling is supported by the evidence 
of doubly peaked electric fields obtained directly 
from beam test measurements and presented in [3]. 
The dependence of the modeled trap concentra- 
tions upon fluence was presented in [4] . In this pa- 
per, we summarize the previous results and inves- 
tigate the temperature dependence of the model. 

This paper is organized as follows: Section 2 de- 
scribes the experimental setup. Section 3 describes 
the carrier transport simulation used to interpret 
the data. The tuning of the double junction model 
and its resulting predictions are discussed in Sec- 
tion 4. The temperature dependence of the data 
and model are summarized in Section 5. The con- 
clusions are given in Section 6. 



2. Experimental setup 

The measurements were performed in the H2 
beam line of the CERN SPS in 2003/04 using 150- 
225 GeV pions. The beam test apparatus is de- 
scribed in [5]. A silicon beam telescope [6] consisted 
of four modules each containing two 300 |J.m thick 
single-sided silicon detectors with a strip pitch of 
25 |Am and readout pitch of 50 |a.m. The two detec- 
tors in each module were oriented to measure hori- 
zontal and vertical impact coordinates. A pixel hy- 
brid detector was mounted between the second and 
third telescope modules on a cooled rotating stage. 
A trigger signal was generated by a silicon PIN 
diode. The analog signals from all detectors were 
digitized in a VME-based readout system by two 
CAEN (V550) ADCs and one custom-built flash 



ADC. The entire assembly was located in an open- 
geometry 3T Helmholtz magnet that produced a 
magnetic field either parallel or orthogonal to the 
beam. The temperature of the tested sensors was 
controlled with a Peltier cooler that was capable 
of operating down to -30° C. The telescope infor- 
mation was used to reconstruct the trajectories of 
individual beam particles and to achieve a precise 
determination of the particle hit position in the 
pixel detector. The resulting intrinsic resolution of 
tlic^ beam telescope was about 1 )im. 

The prototype pixel sensors are so-called "n-in- 
n" devices: they are designed to collect charge from 
n"*" structures implanted into n bulk silicon using 
p-spray isolation. All test devices were 22x32 ar- 
rays of 125x125 /Ltm^ pixels that were fabricated 
by CiS. The substrate, produced by Wacker, was 
285 /im thick, n-doped, diffusively-oxygenated 
float zone silicon of orientation (111), resistivity 
3.7 kfi-cm and oxygen concentration in the or- 
der of 10^^ cm~^. Individual sensors were diced 
from fully processed wafers after the deposition of 
under-bump metalization and indium bumps. A 
number of sensors were irradiated at the CERN 
PS with 24 GeV protons. The irradiation was 
performed without cooling or bias. The delivered 
proton fluences scaled to 1 MeV neutrons by the 
hardness factor 0.62 [7] were 0.5x10^^ Ucq/cm^, 
2x10" Ueq/cm^ and 5.9x10" u^Jcrn^. Afl sam- 
ples were annealed for three days at 30° C. In order 
to avoid reverse annealing, the sensors were stored 
at -20°C after irradiation and kept at room tem- 
perature only for transport and bump bonding. 
All sensors were bump bonded to PSI30/AC30 
readout chips [8] which allow analog readout of 
all 704 pixel cells without zero suppression. The 
PSI30 settings were adjusted to provide a linear 
response to input signals ranging from zero to 
more than 30,000 electrons. 



3. Sensor simulation 

The interpretation of the test beam data relies 
upon a detailed sensor simulation that includes the 
modeling of irradiation effects in silicon. The sim- 
ulation, PIXELAV [2,9,10], incorporates the foUow- 
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ing elements: an accurate model of charge deposi- 
tion by primary hadronic tracks (in particular to 
model delta rays); a realistic 3-D intra-pixel elec- 
tric field map; an established model of charge drift 
physics including mobilities, Hall Effect, and 3-D 
diffusion; a simulation of charge trapping and the 
signal induced from trapped charge; and a simu- 
lation of electronic noise, response, and threshold 
effects. The intra-pixel electric field map was gen- 
erated using TCAD 9.0 [11] to simultaneously solve 
Poisson's Equation, the carrier continuity equa- 
tions, and various charge transport models. A fi- 
nal simulation step reformatted the simulated data 
into test beam format so that it could be processed 
by the test beam analysis software. 

The effect of irradiation was implemented in the 
TCAD simulation by including two defect levels in 
the forbidden silicon bandgap with opposite charge 
states and trapping of charge carriers. The model, 
similar to one proposed in [12], is based on the 
Shockley- Read-Hall (SRH) statistics and produces 
an effective space charge density peS from the trap- 
ping of free carriers in the leakage current. The ef- 
fective charge density is related to the occupancies 
and densities of traps as follows, 



PeflF = e [Nnfo - Na/a] + Pdopants 



(1) 



where: Nn and Na are the densities of donor and 
acceptor trapping states, respectively; fu and fA 
are the occupied fractions of the donor and accep- 
tor states, respectively, and Pdopants is the charge 
density due to ionized dopants (describes the re- 
sistivity of the material before irradiation). The 
donor and acceptor occupancies are related to the 
trap parameters by standard SRH expressions 



fD = 



fA = 



Vh 



(2) 



where: v^. and Vh are the thermal speeds of elec- 
trons and holes, respectively; erf, are the elec- 
tron and hole capture cross sections for the donor 
trap; cr^, are the electron and hole capture cross 
sections for the acceptor trap; n, p are the densities 
of free electrons and holes, respectively; m is the 



intrinsic density of carriers; E'd, Ea are the acti- 
vation energies (relative to the mid-gap energy) of 
the donor and acceptor states, respectively. Note 
that the single donor and acceptor states model the 
effects of many physical donor and acceptor states 
making the two-trap model an "effective theory" . 

The physics of the model is illustrated in Fig. 1. 
The space charge density and electric field are plot- 
ted as functions of depth z for a model tuned to re- 
produce the $ = 5.9 X lO^^UeqCm"^ charge collec- 
tion data at 150V bias. The SRH process produces 
electron-hole pairs more or loss uniformly across 
the thickness of the sensor. As the electrons drift 
to the n+ implant, the total electron current in- 
creases as z decreases. The hole current similarly 
increases with increasing z. Trapping of the mobile 
carriers produces a net positive space charge den- 
sity near the p"*" backplane and a net negative space 
charge density near the n+ implant. Since positive 
space charge density corresponds to n-type doping 
and negative space charge corresponds to p-type 
doping, there are p-n junctions at both sides of 
the detector. The electric field in the sensor follows 
from a simultaneous solution of Poisson's equa- 
tion and the continuity equations. The resulting 
2;-component of the electric field varies with an ap- 
proximately quadratic dependence upon z having 
a minimum at the zero of the space charge density 
and maxima at both implants. A more detailed de- 
scription of the double junction model and its im- 
plementation can be found in [2]. 
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Fig. 1. The space charge density (sohd line) and electric 
field (dashed line) at T = — 10°C as functions of depth in 
a two-trap double junction model tuned to reproduce the 
$ = 5.9x lO-^^neqcm"^ charge collection data at 150V bias. 
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4. Model tuning and results 

Charge collection across the sensor bulk was 
measured using the "grazing angle technique" [15]. 
As is shown in Fig. 2, the surface of the test sen- 
sor is oriented by a small angle (15°) with respect 
to the pion beam. Several samples of data were 
collected with zero magnetic field and at temper- 
ature of — 10°C and — 25°C. The charge measured 
by each pixel along the y direction samples a dif- 
ferent depth z in the sensor. Precise entry point 
information from the beam telescope is used to 
produce finely binned charge collection profiles. 



sity of ionized acceptor impurities cannot repro- 
duce the measured profiles. 
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Fig. 2. The grazing angle technique for determining charge 
collection profiles. The charge measured by each pixel along 
the y direction samples a different depth z in the sensor. 



The charge collection profiles for a sensor irradi- 
ated to a fluence of $ = 5.9 x 10^** Ucq/cm^ and op- 
erated at a temperature of — 10°C and bias voltages 
of 150V and 300V are presented in Fig 3. The mea- 
sured profiles are shown as solid dots and the simu- 
lated profiles are shown as histograms. In order to 
investigate the applicability of the traditional pic- 
ture of type-inverted silicon after irradiation, the 
simulated profiles were generated with electric field 
maps corresponding to two different effective den- 
sities of acceptor impurities. The full histograms 
are the simulated profile for A'cff = 4.5x10^^ cm"-^. 
Note that the 300V simulation reasonably agrees 
with the measured profile but the 150V simulation 
is far too broad. The dashed histograms show the 
result of increasing A^off to 24 x 10^^ cm^"^. At this 
effective doping density, the width of the simulated 
peak in the 150V distribution is close to correct 
but it does not reproduce the "tail" observed in 
the data at large y. The 300V simulated distribu- 
tion is far too narrow and the predicted charge is 
lower than the data (note that the profiles are ab- 
solutely normalized). It is clear that a simulation 
based upon the standard picture of a constant den- 



Fig. 3. The measured and simulated charge collection pro- 
files for a sensor at T = — 10°C irradiated to a fluence of 
$ = 5.9 X 10^"* neq/cm^. The profiles measured at bias 
voltages of 150V and 300V are shown as solid dots. The 
full histograms are the simulated profiles for a constant 
effective doping A^eft = 4.5 X lO^'^ cm~^ of acceptor im- 
purities. The dashed histograms are the simulated profiles 
for a constant effective doping A^^ff = 24 X 10^^ cm~^. 

The same measured profiles and those from bias 
voltages of 200V and 450V are shown in Fig. 4. 
They are compared with simulations based upon 
the electric field produced by the two trap model. 
The model has six free parameters (A^, A^^, cr^, 
cr,f , <jf, cr^) that can be adjusted. The activation 
energies are kept fixed to the values of [12]: Ejj — 
Ev + 0.48 eV, Ea = Ec - 0.525 eV where Ey and 
Ec are the energies of the valence and conduction 
band edges. The electric field map produced by 
each TCAD run is input into pixelav. The electron 
and hole trapping rates, Fe and F^, are also inputs 
to PIXELAV and are treated as constrained parame- 
ters. Although they have been measured [13], they 
are uncertain at the 20% level due to the fiuence 
uncertainty and possible annealing of the sensors. 
They are therefore allowed to vary by as much as 
±20% from their nominal values. The donor con- 
centration of the starting material is set to 1.2 x 
10^^ cm~'^ corresponding to a full depletion volt- 
age of about 70 V for an unirradiated device. Be- 
cause each model iteration took approximately two 
days, it was not possible to use standard statisti- 
cal fitting techniques. The parameters of the dou- 
ble junction model were systematically varied and 
the agreement between measured and simulated 
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charge collection profiles was judged subjectively. 
The "best fits" shown in this paper are probably 
not true likelihood minima and the calculation of 
eight parameter error matrices is beyond available 
computational resources. Adequate agreement was 
achieved by setting the ratio of the common hole 
and electron cross sections Ghl<Je to 0.25 and the 
ratio of the acceptor and donor densities Na/Nd to 
0.40. There is a range of parameters in the N^-ae 
space that produces reasonable agreement with the 
measured profiles. The range is shown in Fig. 5a as 
the solid line in the logarithmic space. If the donor 
density becomes too small {Nu < 20 X 10^'* cm~^), 
the 150V simulation produces too much signal at 
large z. If the donor density becomes too large 
{Nd > 50 X 10" cm-3), the 300V simulation pro- 
duces insufficient signal at large z. Since the simu- 
lated leakage current varies as /loak aeNn, differ- 
ent points on the allowed solid contour correspond 
to different leakage current. Contours of constant 
leakage current are shown as dashed curves and 
are labeled in terms of the corresponding damage 
parameter a where ao = 4 x 10"^'' A/cm is the 
expected leakage current [14]. It is clear that the 
simulation can accommodate the expected leakage 
current which is smaller than the measured current 
by a factor of three. The same choice of parame- 
ters can also account for the observed rate of signal 
trapping [2]. 

The simulation describes the measured charge 
collection profiles well both in shape and normal- 
ization. The "wiggle" observed at low bias voltages 
is a signature of the doubly peaked electric field 
shown in Fig. 1. The relative signal minimum near 
y = 700 yum (see Fig. 4) corresponds to the mini- 
mum of the electric field ^;-component, Ez, where 
both electrons and holes travel only short distances 
before trapping. This small separation induces only 
a small signal on the n+ side of the detector. At 
larger values of y, increases causing the elec- 
trons drift back into the minimum where they are 
likely to be trapped. However, the holes drift into 
the higher field region near the p+ implant and 
are more likely to be collected. The net induced 
signal on the n+ side of the detector therefore in- 
creases and creates the local maximum seen near 
y = 900 /xm. 

The charge collection profiles at T = — 10°C 
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Fig. 4. The measured charge collection profiles at a tem- 
perature of -10°C and bias voltages of 150V, 200V, 
300V, and 450V are shown as solid dots for a fluence of 
5.9 X 10^* Ueq/cm^. The two-trap double junction simula- 
tion is shown as the solid histogram in each plot. 
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Fig. 5. The allowed region in the Njj-ae space for the 
best fit 5.9 X lO"'^ iicq/cm'^ model is shown as the solid 
line. Contours of constant leakage current are shown as 
dashed curves and arc labeled in terms of the corresponding 
damage parameter a where ao = 4 x 10"'^'^ A/cm is the 
expected leaJjage current [14]. 

for sensors irradiated to fluences of $ = 0.5 x 



10 



14 



j/cm^ and <[> = 2 x 10 Ueq/cm^ and op- 



erated at several bias voltages are presented in 
Fig. 6(a-c) and Fig. 6(d-g), respectively. The mea- 
sured profiles, shown as solid dots, are compared to 
the simulated profiles, shown as histograms. Note 
that the "wiggle" is present at low bias even at 
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$ = 0.5 X 10^* neq/cm^ which is just above the 

"type-inversion" fluence. This suggests that a dou- 
bly peaked field is present even at rather small flu- 
ences. 
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Fig. 6. Measured (full dots) and simulated (histogram) 
charge collection profiles for sensors irradiated to fiuences of 
* = 0.5 X 10" neq/cm2 (a-c) and of * = 2 x lO" Ueq/cm^ 
(d-g), at T = — 10°C and several bias voltages. 

The double junction model can provide a rea- 
sonable description of the lower fluence charge col- 
lection profiles using the parameters obtained with 
the fitting procedure shown in Table 1. We ob- 
serve that the donor trap concentration increases 
more rapidly with fluence than does the accep- 
tor trap concentration. The ratio between acceptor 
and donor trap concentrations is 0.76 at the lowest 
fluence and decreases to 0.40 at 5.9x10^^ Ueq/cm^. 
In addition, the fits exclude a linear dependence 
of the trap concentrations with the irradiation flu- 
ence. At $ = 5.9x10^^ neq/cm^ the cross section 
ratio ah/de is sot to 0.25 for both donor and accep- 
tor traps while at lower fluences we find cr^/cr^ = 
0.25 and c^/CTf = 1 for the acceptor and donor 
traps, respectively. The simulated leakage current 
is approximately linear in fluence, but the ratio 
Na/Nd is clearly not constant. This may be a con- 
sequence of the qtiadratic fluence scaling of one or 
more di-vacancy states or it may reflect the fact 
that the two trap model with the particular choice 
of activation energies docs not accurately model 
the dependence of the trap occupancies on leak- 
age current. The allowed Njj-cre parameter spaces 
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Fig. 7. The z-component of the simulated electric field at 
T = — 10°C resulting from the model best fit is shown 
as a function of z for a sensor irradiated to fluences of 
* = 0.5 X 10" Ueq/cm^ (a) and * = 2 x lO" Ueq/cm^ (b). 



for the lower fluence models arc much more con- 
strained than in the $=5.9 x 10^** neq/cm^ case anc 
predict the expected leakage current. 



* [10" neqCm-2] 0.50 2.0 5.9 
Na [10" cm-3] 1.9 6.8 16 



Nr, [10" cm- 



2.5 10 40 



at^'^ [10-15 cm^] 6.60 6.60 6.60 



[10-15 



1.65 1.65 1.65 
6.60 6.60 1.65 



Table 1 

Double trap model parameters extracted from the fit to 
the data. 

The z-component of the simulated electric field, 
Ez, is plotted as a function of z in Fig. 1 for $ = 
5.9 x lO^'^ noq/cm^ and in Fig. 7 for $ = 0.5 x 
10" Ueq/cm^ and $ = 2 X 10" Ugq/cm^. At $ = 
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Fig. 8. The simulated space charge density at T = — 10°C 
as a function of the z coordinate for fiuences of 
* = 0.5 X 10" Ueq/cm^ and * = 2 X lO" Ueq/cm^. 
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5.9 X 10^"' ncq/cm^, the field profile has a min- 
imum near the midplane of the detector which 
shifts toward the p+ implant at lower fluences. 
The field has maxima at the detector implants 
as discussed in Section 3. Figure 7(a) shows that 
a doubly peaked electric field is necessary to de- 
scribe the measured charge collection profiles even 
at the lowest measured fluence which is just be- 
yond the "type inversion point" . The dependence 
of the space charge density upon the z coordinate 
is shown in Figures 1 and 8. Before irradiation the 
sensor is characterized by a constant and positive 
space charge density of 1.2 x 10^^ cm"-^ across the 
sensor bulk. After a fluence of 0.5 x 10^'* Ucq/cm^ 
the device shows a negative space charge density of 
about —1 X 10^^ cm~'^ for about 70% of its thick- 
ness, a compensated region corresponding to the 
Ez minimum and a positive space charge density 
close to the backplane. The space charge density 
and electric field near the p+ implant increase with 
fluence. The space charge density is not linear in 
z due to the variation of the carrier drift velocities 
with the electric fields. 



5. Temperature dependence 
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Fig. 9. Measured (full dots) and simulated (histogram) 
charge collection profiles at T = — 25°C and sev- 
eral bias voltages for sensors irradiated to fluences of 
* = 2.0x10^* neq/cm^ (a-d) andof <t> = 5.9xl0l-* Ucq/cm^ 
(e-g). 



The temperature dependence of the charge col- 
lection profiles was studied by accumulating data 
at r = — 25°C. The pixelav simulation includes 
temperature dependent mobilities, diffusion, and 
trapping rates. The tcad calculation of the elec- 
tric field map is also based upon temperature de- 
pendent quantities including the bandgap energy 
and SRH lifetimes. The T = -25°C charge collec- 
tion profiles for the $ — 2.0 x lO^"' Ucq/cm'^ and 
$ = 5.9 X 10""^* Ucq/cm^ sensors are compared with 
the simulation in Fig. 9. It is clear that the simu- 
lation correctly tracks the temperature-dependent 
variations in the measured profiles. 

The effect of temperature on the z-component 
of the simulated electric field at $ = 5.9 x 
10"'^^ Uoq/cm^ is shown in Fig. 10 for bias voltages 
of 150V and 300V. It is clear that decreasing the 
temperature also decreases the fields on the p-|- 
side of the sensor and increases them on the n+ 
side. 
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Fig. 10. The simulated 2-component of the electric field as a 
function of the z coordinate at the <J> = 5.9 x 10^^ neq/cm^ 
fiuence for temperatures T = -10°C and T = -25°C. The 
field profiles are shown for bias voltages of 150V and 300V. 
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6. Conclusions 



In this paper we have shown that doubly peaked 
electric fields are necessary to describe grazing- 
angle charge collection measurements of irradiated 
silicon pixel sensors. A model of irradiated silicon 
based upon two defect levels with opposite charge 
states and the trapping of charge carriers can be 
tuned to produce a good description of the mea- 
sured charge collection profiles in the fluence range 
from 0.5x10" Uoq/cm^ to 5.9x10" n^Jcm^. The 
model correctly predicts the variation in the pro- 
files as the temperature is changed from — 10°C to 
-25°C. 

The doubly peaked electric field profiles have 
maxima near the implants and minima near the 
detector midplane. This corresponds to negative 
space charge density near the n"*" implant and and 
positive space charge density near the p"*" back- 
plane. We find that it is necessary to decrease the 
ratio of acceptor concentration to donor concentra- 
tion as the fluence increases. This causes the elec- 
tric field profile to become more symmetric as the 
fluence increases. The effect of decreasing the tem- 
perature has the opposite effect of suppressing the 
fields on the p-l- side of the sensor and increasing 
them on the n-|- side. 

The measured charge collection profiles of irra- 
diated sensors are inconsistent with the linearly- 
varying electric fields predicted by the usual de- 
scription based upon a uniform effective doping 
density. This suggests that the correctness and the 
physical significance of effective doping densities 
determined from capacitance-voltage measure- 
ments are quite unclear. In addition, we remark 
that the notion of partly depleted silicon sensors 
after irradiation is inconsistent with the measured 
charge collection profiles and with the observed 
doubly peaked electric fields. 

The charge-sharing behavior and resolution 
functions of many detectors are sensitive to the 
details of the internal electric field. A known 
response function is a key element of any recon- 
struction procedure. A working effective model 
will permit the detailed response of these detec- 
tors to be tracked as they are irradiated in the 
next generation of accelerators. 
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